Parity nonconservation in deuteron photoreactions 
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We calculate the asymmetries in parity nonconserving deuteron photodisintegration due to cir- 
cularly polarized photons yd — > np with the photon laboratory energy ranging from the threshold 
up to 10 MeV and the radiative capture of thermal polarized neutrons by protons np — > -yd. We 
use the leading order electromagnetic Hamiltonian neglecting the smaller nuclear exchange currents. 
Comparative calculations are done by using the Reid93 and Argonne wis potentials for the strong 
interaction and the DDH and FCDH "best" values for the weak couplings in a weak one-meson 
exchange potential. A weak JVA transition potential is used to incorporate also the A(1232)-isobar 
excitation in the coupled-channels formalism. 



INTRODUCTION 



Due to the incessant presence of the strong interaction, the hadronic weak interaction in the quark flavour conserving 
sector of the Standard Model is not completely understood. Even though the weak interaction between quarks is well 
qq " known at high energies, its properties are hard to extract in the nonperturbative regime of quantum chromodynamics 
(QCD). This is due to the complicated structure of nuclear systems along with coinciding dynamics of QCD. It 
seems an almost hopeless challenge to distinguish directly about the size of seven orders of magnitude smaller weak 
interaction effects from those of QCD. Fortunately, the weak interaction leaves a unique signature in the form of 
parity nonconservation (PNC), which provides a tiny but non- vanishing observable. 

For the past three decades the PNC calculations between nucleons have been based on the use of the DDH 1] 
potential. Unfortunately, the potential requires the knowledge of several weak meson-nucleon coupling constants, 
which are still to date inadequately known. Today's contemporary attempts to determine PNC amplitudes more 
and more often harness QCD based chiral perturbation theory (%PT). Unlike the phenomenological meson-exchange 
model, xPT provides a systematic and model-independent way to study hadronic reactions at low energies. The 
updated xPT based approach leads to a systematic expansion of PNC amplitudes in terms of low-energy constants 
(LECs) which have a straightforward correspondence with the weak DDH meson-nucleon couplings, see e.g. refs. 

Out of all the possible candidates for asymmetry and polarization observables, in this paper we focus on the 
asymmetries associated with PNC deuteron photodisintegration by circularly polarized photons jd — > np and thermal 
polarized neutron capture np — > "fd. It is further to be noted that the photon polarization of np — ¥ yd at threshold 
equals the photon asymmetry of the time-reversed reaction jd — > np at threshold. The threshold behaviour of 
the asymmetry/polarization in reactions 7c? O np can be shown to be insensitive to the tt- meson exchange, which 
represents the long-range AI — 1 part of the PNC interaction. Therefore, the threshold region is essentially dominated 
by the exchanges of heavy mesons (p, u) and thus also the relatively long-ranged A-excitation could occur more 
clearly highlighted than what it would if it appeared in a background where the pion is more intensely present. To our 
knowledge this effect has only been checked in the form of exchange currents in refs.[l,Q. Contrary to the reactions 
yd -H- np at the threshold, the PNC 7r-exchange is predominant in the low-energy reaction np — > "/d. 

There are various theoretical works on the PNC reactions yd —> np [3, 134l3|. np — > -yd [l2Hl9l |. and np —> 7c? 
[1, IE EH EH EH EH H3-HH]. The calculations are typically carried out exploiting the old meson-exchange picture, 
aside from some of the recent works, which apply the modern state-of-art techniques, e.g. such as the pionless 
effective field theory EFT(x) and heavy-baryon chiral perturbation theory HBxPT. The results more or less agree in 
the threshold region. The energy regime of several MeV above the threshold of the deuteronphotodisintegration is 
investigated in refs. pi [7rlTo|. Again, in that regime, the results are similar except in rcf. 7], which differs by an 
exceptionally large pion contribution. The difference is discussed in detail in ref. [1(|. Up to date, for the reaction 
•yd —> np, there exist only two experimental data points: (7.7 ± 5.3) x 10~ 6 and (2.7 ± 2.8) x 1CU 6 at the photon 
energies of 3.2 and 4.1 MeV respectively [29| . The latest photon polarization measurement for the inverse reaction 
np — !• t*g? gives the value of (1.8 ± 1.8) x 10 -7 [3(|. The data from the 1980's for both the reactions are consistent 



'Electronic address: tero.partanen@helsinki.fi 
tElectronic address: jouni.niskancn@hclsinki.fi 



2 



with zero with rather a poor precision and, therefore, new and more accurate experimental data, e.g. [31L l32j. would 
be appreciated. 

The asymmetry of the reaction rip — > jd was also measured previously in the 70's and 80's, resulting in the 
insufficiently accurate values of (0.6 ± 2.1) x 10 -7 [33[ and —(1.5 ± 4.8) x 10 -8 (34J. However, there is currently 
an ongoing experiment (NPDGamma) on the reaction rip — > ~fd with cold neutrons at LANSCE and SNS with a 
preliminary result, which sets the asymmetry within the limits of (—1-2 ± 2.1(stat.) ± 0.1(sys.)) x 10~ 7 [35]. The 
oncoming measurements of the NPDGamma experiment aim to improve the accuracy up to a level of 20% of the 
typical theoretical prediction —5 x 10 -8 , which employs the DDH "best value" for the weak irNN coupling constant 
h^k \ The asymmetry is straightforwardly proportional to the and therefore the current experiment is about to 
shed some light on the uncertain value of the coupling. 

In this paper we study PNC in the above reactions using two modern phenomenological strong potentials and 
post-DDH weak couplings. Furthermore, we want to estimate the size of the aforementioned A effect, which was 
found to be significant in PNC elastic pp scattering at higher energies [36l . 113] • As in these works our calculation is 
carried out within the framework of the coupled channels meson-exchange model and hence, according to the common 
practice, we utilize the DDH potential [l[ as the starting point for the PNC it-, uj-, and p-exchanges extending to use 



the weak couplings of ref. [38[ consistent with the presence of the A. We take account of exchange currents only in 
the extent they exist when the Siegert's theorem [39j is applied. They have been considered more explicitly e.g. in 
refs. [H. IH Hoj|. In practise we use the electromagnetic Hamiltonian under the dipole approximation, which allows the 
PNC deuteron breakup to have four pn continuum channels 1 So, 3 Po, 3 Pi, and 3 p2 — 3 P2 in the expansion up to P 
waves. 

This paper is organized as follows. In sect. UH we give the appropriate forms of the electromagnetic Hamiltonian, 
the scattering and deuteron wave functions, and the PNC one-meson exchange JVJV and TVA-transition potentials. 
The forms of the spin observables A^ and An respectively of the reactions 7c? — > np and rip — > jd are also given. 
Summary of the results and conclusions are presented in sect. IIII1 

II. THEORY 
A. Interactions 

The observable asymmetries in the photoreactions arise from an interference between the strong and weak interac- 
tions giving rise to simultaneous presence of the photomagnetic and photoelectric effects. The interference appears 
when these two parallel processes share the same final state quantum numbers, and the PNC observables are obtained 
in terms of products of PC and PNC partial wave amplitudes. This is illustrated in fig. [TJ where the mechanisms 
(TV TV, iVA with PC and PNC forces) are shown alongside with the quantum numbers of the possible states up to 
P-wave final nucleons. The electromagnetic perturbing Hamiltonian H e . m . = H^ N + ff^f + ff™ which takes care 
of the disintegration and formation of the deuteron is considered in the dipole approximation e ±lk i- r « 1. 

The Hamiltonian for absorption of a photon may be written as 



H c \u. = -iJ^e^x ■ + + ^(C 1 + W + (P> + *>n)°i + W(T?Si + h.c.)) x fe 7 ] , (1) 

i 

where e^x (A = ±1) is the circular polarization vector of the incoming photon, oj 1 the center-of-mass energy of the 
photon, /i s = 0.88 and fi v = 4.71 the isoscalar and isovector magnetic moments of nucleons, a = e 2 the fine-structure 
constant, and M = 939 MeV the average nucleon mass. The 77VA-vertex can occur in the presence of the Ml and E2 
transitions. We consider only the dominant Ml multipole and neglect the small E2 effect. The value of the transition 
magnetic moment is given by the quark model as fj,* = f*/i v /2f = 3\Z2/j, v /5 |40l |. The nucleon-Delta spin and isospin 
transition operators are denoted as S and T z [4l|. 

In the presence of possible channel coupling the final state scattering wavefunctions are of the form 

& s (k,r) = ^EE i L {LM L SMs\JM)Yl ML {k)U^7\k,r)yyS\f)\T%, (2) 

T k'JMLM l 



where the superscript " (— )" on the wavefunctions refers to the incoming wave boundary conditions. The |T0) are the 

\) LS ( 



relevant isospin states for the pn-interaction, 3^jm(^) the eigenfunctions of the coupled total angular momentum, and 



' '{k, r) the complex- valued radial wavefunctions. The quantum numbers LST (relative orbital angular momentum 
L, total spin S, and total isospin T), which may be changed by the nuclear forces, are denoted for brevity by k which 
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FIG. 1: Graphical representation for the np and the direct iVA contributions in the deuteron photodisintegration. The dashed 
line represents the strong interaction, the dotted line the weak interaction, and the wavy line the incoming photon. For brevity 
the tensor coupled scattering states 3 5*i — 3 Di and 3 P 2 — 3 7 7 2 are denoted by 3 Ci and 3 C 2 respectively. Similarly, the d stands 
for the standard PC deuteron 3 Si — 3 Di state. 



also refers to the main final wavefunction, whereas k! labels the "small" component generated from it. The J and M 
relate to the total angular momentum and they are good quantum numbers under the nuclear forces. The calculated 
matrix elements '~* ) (fe; SMsT\H^ m \Md)v will be expressed in terms of the wavefunctions with the outgoing boundary 
conditions "(+)", i.e. U^t\k,r) =U J J~ ] * '(k,r). 

The deuteron wavefunction consists of a superposition of six relevant components: the standard PC NN ( 3 <Si, 3 Di) 
and the tiny PNC NN ( 1 P 1 , 3 Pi) and NA ( 3 P 1 , 5 P 1 ) states and may be written as 

<U ( r ) = E ^T^yfa 3 / ( f )\T d 0) , (3) 

with the normalization dr^i \^i\ 2 = 1- The PC bound wavefunctions are real- valued and their PNC partners 
imaginary- valued. Both the scattering and bound state wavefunctions are obtained by solving coupled Schrodinger 
equations. The radial Schrodinger equation for the parity-admixed interaction reads 

= 2„yj / Jfi r y™t(r)(T"0jl> po (r) + l>™ o (r)|T'0)^OTW„ J «(*,r), (4) 

k" 

where for V PC (r) the digonal NN potential is taken as the phenomenological updated Reid soft core (Reid93) [42[ 
and Argonne v\% (Avis) [HI potentials and respectively the AN as the mass difference Ma — M to generate the 
strong correlations. The other potentials are defined in eqs. (0 and (l8|)- (fTTj) . In the presence of AN channels the 
TV TV potentials must be modified to avoid doubly counting the attraction from the channel coupling guaranteeing the 
phase equivalence. It is crucial for the result that it is done with high precision. The reduced masses fi for the initial 
NN and AiV states are respectively fi NN = M/2 and pi NA = MM a/ (M + M A ). 
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TABLE I: The weak aNN couplings ha and aNA couplings h a . The set of the first three weak couplings are the DDH 
"best values" [lj and the following five the FCDH "best values" [38l ]. The weak couplings are given in units of 10~ 7 . 
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For the weak NN interaction we use the DDH weak one-meson exchange nucleon-nucleon potential ref. [T[ . In the 
case of pn-system, the relevant part of the potential is 
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where C Q = [— iV,Y Q (r)] and 0+ = {— iV,Y Q (r)}, with cv = 7T, p, a; are respectively the commutator and anticom- 
mutator in which the radial functions are 
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if form factors are not used. In case monopole form factors of the type (A 2 — rn 2 i )/(q 2 + m 2 a ), which we use here, are 
included in vertices, the modified Yukawa functions take the form 
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(7) 



Anr An \r 2A a 

In eq. ([5]) we have neglected the term proportional to hp because of the smallness and vagueness of the coupling 
and also the irrelevant oc (t* +r 2 ) terms, which do not contribute in pn interaction. 

The PNC transition potential (NN O AN) may be derived from the vertex interaction Hamiltonians (|A1[) - (|A8[) of 
Appendix El resulting in the potentials (^-(fTO)) for p-,cv-, and 7r-exchanges respectively: 
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FIG. 2: Diagrammatic representation of the isospin change AI in the process of the PNC reaction yd — > np through the direct 
A-channel. The subscript denotes the isospin state. 
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where g* a = y/72/25g a ( a = Tr,p) are the strong meson- TV A couplings following from the quark model. It may be 
noted that in the case of pion-exchange the A is generated only at the strong vertex. The ir and p -mediated strong 
transition potential of the standard form [41| is 



KHr) = gf§Ti • r 2 (5r • V)(<x 2 • V)Y.(r) 



4M 2 



(1 + XpfTi ■ r 2 (S 1 x V)(<7 2 x V)y p (r) + (1 o 2) + h.c, 



(11) 



involving a spin-spin and tensor part. Note that the strong NA transition potential (|11[) is more singular than r~ 2 and 
thus necessarily requires regularization. Therefore, in order to be thoroughly consistent, we always use the modified 
Yukawa functions ([7]) in the presence of the A. 

The whole process of the PNC photodisintegration yd — > np is compelled to change the initial (PC bound) isosinglet 
state to the final (continuum) isovector state. The disintegration isospin transition operator allows only the AI — 1 
transitions, with (10|Ti Z |00) = (20|Ti 2 1 10) = \/2/3 being the only nonzero matrix elements, where £ = 1,2 labels the 
particle. Only the isovector mesons can couple to the PC N A-vertex, which automatically excludes the w-exchange 
in such amplitudes. The weak 7r-exchange is also excluded, since the total isospin would not change in processes via 
A channels. For the same reason there is no contribution coming from the PNC A^Aa; vertex. In general, all the 
structures of the PNC transition potential related to the isospin AI = 1 change are zero, and thus only the p-exchange 
occurs in the presence of the A-channel, see fig. [21 

The potential (jHJ has only one relevant isospin transition amplitude lolin scattering and basically two 0f}0 and 
O 2 in the deuteron, from which only the latter one is nonzero. Following the usual practice in similar calculations, 
we also choose to use the DDH couplings in order to be comparable with the corresponding works. Unfortunately, 
there are no published weak aNA (a = ir,p,Lu) couplings corresponding to the weak DDH aNN couplings. The 
" best" values of the weak aTVA-couplings have been evaluated in newer analyses [38, H3| ■ To extract the pure A effect 
we shall use the FCDH values [38| for the needed weak aA^A-couplings, even though they would not necessarily be 
entirely consistent with the other DDH couplings. However, we also study the effect of using the couplings from the 
consistent analysis of ref. (38j . 

It may be noted that, as a consequence of coupled-channel dynamics, in addition to the one-meson exchanges 
depicted in fig. [1] the A channels may also include possible higher order corrections, which are naturally taken into 
account as correlation effects. For instance, at the threshold, the leading A contributions originate from the once- 
iterated meson exchange diagrams presented in fig. [3l where eq. (|llj) is employed to take care of the strong iVA 
transitions. 



B. Observables 



As a first test we calculate the total photon absorption cross-section a 1 = <j+ + <r_ for the reaction yd — > np, 
obtained by Fermi's golden rule summing over the photon helicities and the final spin states and averaging over the 
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FIG. 3: Second order A corrections at the threshold. 



two polarization directions of the photon and the three possible spin projections of the deuteron. Due to the identity 
of the final state particles (in the isospin formalism) the total cross section is also divided by two. Thus we have 



• W 



(\£^\ 2 + \Mfj 2 



(12) 



where M.^. and £*jf are respectively the general reduced matrix elements of the magnetic and electric transitions 

given in Appendix [BJ The isovector -M 3S transition is the only non-negligible Ml contribution at low energies and 
is furthermore non- vanishing only in the threshold domain. The result is in good agreement with experimental data, 
as seen in figs. H]and[SJ In the latter figure we have used the reciprocity relation 



= 2 V ~fc 



'7! 



(13) 



for the inverse reaction np —> jd. 

The relevant asymmetry observable for the reaction 7c? — ¥ np may be expressed in terms of the deuteron photodis- 
intcgration hclicity cross sections o\ as 



E KKd j(\MZ J d \> + \£g\* + \£^\ 2 + \M^ 



(14) 



Here and later on, the parity admixed (weak) wavefunctions and amplitudes generated by the weak nuclear force will 
be tilded for clarity. Naturally, the absolute squares of the weak amplitudes in the denominator of eq. (| 14[) may as 
well be ignored due to their diminutive size. Specializing to low partial waves they reduce to the results given e.g. in 
refs. [10( and [151 ]. The asymmetry may also be written following the notation fixed in fig. [T]as 
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FIG. 4: The photodisintegration cross-sections of a deuteron for El, Ml, and their sum given by eq. (|12[) . The Reid93 potential 
is employed in the left-hand side figure and the Avis potential in the right-hand side figure. The data points (□) for Ml are 
taken from ref. (4f| and (+) for the total cross-section from ref. [4(J. The w 7 is the incident photon laboratory energy. 
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FIG. 5: The neutron capture cross-sections on a proton given by eq. (|13|) (the Reid93 potential is used in the left and the 
A«i8 right). The data is from refs. [47l. |48|]. The lack of the meson-exchange currents (mainly the one-pion-exchange current 
correction) weakens the Ml of about 10% at low energies [ZjJ. The neutron speed v n (in the units of meters/nanosecond) and 
kinetic energy T n are in the laboratory frame. 



where the index i denotes the final continuum channel. The total asymmetry A^ with S and P waves is achieved 
when k — 4. 

The magnetic dipole effect is dominant at the threshold where the asymmetry (|15|) reduces to a simple form 
Af T -(k -> 0) ~> 2Re[£ 1 /M 1 ], which is explicitly given by eqs. (|B1[) and (|B2[) as 
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(16) 



It is important to note that eq. (fT6| arises from the spin changing PNC interaction, see fig. [TJ and thus does not 
include the PNC pion exchange. Therefore, contributions from heavier vector meson exchanges (and possibly A) are 
maximized. The low-energy limit (1161) coincides with the photon polarization in the time-reversed reaction np — > jd 
for thermal neutrons. 

The asymmetry from the radiative capture rip — > 7c? of longitudinally polarized thermal neutrons in hydrogen is 
also calculated. The appropriate scattering wavefunctions depending on the spin magnetic quantum numbers (m n 
and Trip) of the neutron and proton are obtained by expanding eq. ([2]). The wavefunctions (assuming the z-axis to 
be along the direction of fc) become 
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(17) 



Since this is a time-reversed process to the deuteron disintegration, the amplitudes are of the form 

;m n m p )r'. The asymmetry observable, given in terms of spin differential neutron capture cross sec- 
tions da mrl I dQ with the neutron polarization m n and further in terms of the reduced matrix elements approximated 
for the thermal neutrons, reads An(6) — An cos 8 with 
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(18) 



where 6 is the angle between the momenta of the incident neutron and emitted photon. Figure [6] presents this 
threshold result diagrammatically. Contrary to eq. (|16[) . the An is dominated by the spin conserving weak interaction 
and therefore the PNC one pion exchange prevails. Contrary to 7c? — > np in this reaction the A-corrections are due 
to 7T, p, and to exchanges. 
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FIG. 6: Graphical representation for the PNC reaction np — > 7c! with thermal neutrons. 
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FIG. 7: The PNC deuteron components (divided by the DDH fti ) including the A with the DDH couplings and modified 
Yukawa functions. The Reid93 potential is used in the left panel and the A?;i8 in the right. 



III. RESULTS AND DISCUSSION 



Utilizing the different model complexes outlined in the previous section, we now proceed to calculate the asymmetries 
A^ for the reaction jd — >■ np as a function of photon laboratory energies varying from the deuteron breakup threshold 
to 10 MeV and An in the radiative capture of thermal neutrons (25 meV) np —> -fd. It is also worth noticing that 
the photon polarization of the reaction np — ¥ 7c? with thermal neutrons is in principle the same as the observable 
of the time-reversed reaction at threshold V=j « A^(w^ h '). We also compare in some detail the asymmetries with and 
without the effects of the virtual A-isobar. 

In figs. [7] and [5] we show the PNC deuteron components relevant for the reactions discussed in this paper with 
and without form factors and also including AiV admixture. Two phenomenological potentials are used, the updated 
Reid soft core [42| and Argonne vis [43| potentials. By far, the largest PNC component is 3 Pi, which mainly arises 
from pion exchange as seen in fig. O However, as was seen in sect. [TT] and fig. [I] this does not contribute to PNC 
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FIG. 8: The same as fig. but without the A. The components are compared using both unmodified and modified Yukawa 
functions. The result (unmodified) gained with Avis potential is the same as in ref. (la ]. 
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FIG. 9: Contributions of different mesons in different partial waves without the A using the DDH couplings, the Avis potential, 
and unmodified Yukawa functions. 





0.14 




0.12 




0.1 




0.08 


'e 


0.06 




0.04 


L. 

"O 


0.02 


=3 











-0.02 




-0.04 



i 


i i i 

1p1 _ 




3p1 




A 3p1 - 


_ / 

/.-■-. 

> 


A 5p1 


i i i i 



4 6 
r(fm) 



10 



FIG. 10: The same as fig. [7| The components (divided by the DDH ) with the A«ig potential and the FCDH couplings 



photoabsorption into the lowest 1 Sq partial wave and, consequently, is of minor importance at threshold. The Pi 
component arises from shorter ranged vector meson exchanges, notably p and is smaller by an order of magnitude. 
The pion generated part 3 Pi is relatively model independent (except for the weak pion coupling) whereas in l P\ some 
short-range dependence can be seen, if the PNC potentials are not moderated by a form factor. Also it can be noted 
AN components can be of the same order as 1 P\ although of shorter range. Therefore, it is plausible that their effect 
could be, in principle, appreciable and should be considered. Further, from fig. [5] it is obvious that with the standard 
DDH couplings lo exchange would be negligible and also p in the pion dominated 3 Pi component. Figure [TU] shows 
another choice of weak couplings and the dependence on these is significant. The FCDH parametrization gives smaller 
PNC in comparison with fig. [7J Even though we have compromised to use the FCDH values, the other available weak 
coupling set with A by Desplanques fiij is alike and should give a similar result. Otherwise, if the weak couplings 
are the same, our PNC NN wavefunctions are not in contradiction with those displayed in refs. [1, l9l ITEj]. 

We now consider the asymmetry A^, covering the photon energies from threshold to 10 MeV. Figure ITT1 illustrates 
separately the asymmetry of the "basic model" (by which we mean the commonly used model in the PNC calculations, 
which employs the DDH and strong interaction parameter values of table Q] without form factors) and the effect of the 
form factors and the A employing Reid93 and A^ig potentials. In the case of the "basic model", the asymmetry is 
nearly 30% bigger with Reid93 at threshold. This difference is explained by the distinctly bigger PNC 1 P\ component 
in the deuteron produced by Reid93 potential, as seen in fig. |SJ However, with form factors in the PNC potentials 
the short-range differences are minimized and the results become graphically indistinguishable. 

These points are further featured in tables [II] (last three columns) and IIIII in the threshold limit u>}^ — > 2.22592 
MeV. Short-range correlation differences seem appreciable only in the transition from the PNC deuteron component 
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TABLE II: Contributions of different mesons to the asymmetry A^ = ^2 a A~, in units of 10 8 , at threshold without A-excitation 

using the unmodified Yukawa functions eq. ([6]). The threshold asymmetries are split into two pieces see 

eq. (|16|) . where the indices characterize the El transition: subscripts the initial and superscripts the final states (d for the PC 
part). The Alf = in every case. 



Model 


A' ip oP A 1 Sap AP 


•^■d -^1 p, "^7 


a p i c 
A ° A A~ 


Avis & DDH 
Reid93 & DDH 
A«i 8 & FCDH 
Reid93 & FCDH 


-0.85 3.51 2.66 
-0.84 4.10 3.26 
0.41 1.17 1.58 
0.40 1.36 1.76 


-0.39 0.29 -0.10 
-0.38 0.38 0.00 
-1.00 0.76 -0.24 
-0.99 0.99 0.00 


-1.24 3.80 2.56 
-1.22 4.48 3.26 
-0.59 1.93 1.34 
-0.59 2.35 1.76 



TABLE III: Asymmetries A^ and A^, in units of 10 8 , respectively without and with A-excitation using the modified Yukawa 
functions eq. (0. 



Model 


A 3 Pa a l Sq a 
■^-d p^ ■ n l 


/ p oA A L S A A A 
•^■d -^1 p^ "^7 


Avis & DDH 
Reid93 & DDH 

A«i 8 & FCDH 
Reid93 & FCDH 


-0.69 1.74 1.05 
-0.69 1.83 1.14 
-0.33 0.95 0.62 
-0.33 1.03 0.70 


-0.73 1.69 0.96 
-0.73 1.77 1.04 
-0.27 0.92 0.65 
-0.28 1.00 0.72 



1 P\ to the 1 Sq final state without form factors (table|TT]), while with form factors the difference is hardly 10 %. The 
effect of the form factors themselves in PNC heavy meson exchanges is, however, a dramatic decrease to less than 
half of the original value of the asymmetry A^. In general the FCDH results are much smaller that DDH with and 
without form factors as seen in fig. and Tables. The threshold result = 2.56 x 10~ 8 gained with the "basic 
model" and the A^is potential argrees with most of the existing calculations and perfectly with those of using the 
same model [HI [l5j]. Even though the NA components in the deuteron appeared significant in fig. [3 their effect in 
A^y remains negligible at low energies being about 10% decrease for DDH at threshold and practically null for FCDH. 

Figure Q2] shows the importance of the P-wave continuum states above threshold as a cumulative behaviour of eq. 
(TTSl) . In the same figure are also shown the contributions of different mesons. As the photon energy increases, the 
asymmetry decreases steeply within the energies up to of about 1 MeV above the threshold. This energy region is 
dominated by heavy meson exchanges, mainly the p. Increasing the energy further, the asymmetry quickly settles 
down near to zero continuing its gentle monotonic decline through zero somewhere between 3 and 4 MeV, after which 
the pion starts to take the dominance. 

At low energies the pion contribution to A^ is small and, in our model and expressions, vanishes in the threshold 
limit (there is a negligible contribution higher order in fc 7 arising from the spin current 0, Hoj). This is not so for 




FIG. 11: The total asymmetry given by eq. (|15[) and the DDH couplings with unmodified and modified Yukawa functions 
without (NN) and with the A excitation (NN+NA). The Reid93 left and the A^i 8 right. 
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FIG. 12: The total asymmetry given by eq. I)15[). the Avis potential, and the FCDH couplings with unmodified and modified 
Yukawa functions without (NN) and with the A excitation (NN+NA). 
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FIG. 13: Cumulative behaviour of the asymmetry eq. (|15[) in the left-hand side panel. Single meson tt,uj, and p contributions 
to the asymmetry in the right-hand side panel. Both figures are plotted employing the DDH couplings, unmodified Yukawa 
functions and A«ig potential without A excitation. 



radiative capture of longitudinally polarized neutrons, where it can participate in spin conserving PNC transitions 
3 Si + 3 Di <->• 3 Pi as depicted in fig. [6] and expressed in eq. JTBJ)- As seen in table IIV1 contrary to Ay, now the vector 
meson contributions are negligible. (It may still be of interest to note that now due to the spin conserving couplings 

tJ is more important than p.) PNC in both bound states and continuum is about equal in magnitude and their 
contributions add constructively. Also in this case the value of An = —5.39 x 10~ 8 by the "basic model" and the Avis 
potential is in agreement with results of most authors and especially ref. [i~5| in which the wavefunctions of the same 
model are used. Due to the long range of the pion and low energies the form factor effect is only a 2-3 % decrease in 



TABLE IV: Contributions of different mesons to the An = ^2 a A%, in units of 10 8 , without A-excitation using the unmodified 
Yukawa functions eq. ([6]). Furthermore, as in table[H] the A!% is a sum of PNC scattering Af'p^ and bound Aa^" contributions 
corresponding to eq. (|18[) . 



Model 


A* p i r AZ 

"^3 p 1 "^3 d ^ n 


.dp A^PlP AP 
-^3 Pl 


•^ i 3p 1 * Al 3(^ 1 "^n 


A d ~ /4 dpi A~ 

"^3 p 1 ^3 Q 1 SMl 


Avis & DDH 
Reid93 & DDH 
Awi 8 & FCDH 
Reid93 & FCDH 


-2.85 -2.51 -5.36 
-2.85 -2.37 -5.22 
-1.67 -1.47 -3.14 
-1.68 -1.39 -3.07 


0.00 0.00 0.00 
0.02 0.00 0.02 
0.00 0.00 0.00 
0.01 0.00 0.01 


-0.01 -0.02 -0.03 
-0.03 -0.01 -0.04 
-0.02 -0.05 -0.07 
-0.02 -0.05 -0.07 


-2.86 -2.53 -5.39 
-2.86 -2.40 -5.26 
-1.69 -1.52 -3.21 
-1.69 -1.44 -3.13 
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Aft- The A excitation contributes another 5 % in both continuum and bound states, but these effects cancel off. 

In summary, we have calculated the asymmetries A^ and An in polarized photon absorption and in radiative capture 
of polarized neutrons close to threshold. The results are mutually complementary in the sense that Ay is dominated 
by vector meson exchange, while Ah is pion dominated. The dominances are more than an order of magnitude with 
the minor effect being negligible. Also the A effects are small. In addition, very soon above threshold Ay becomes 
probably too small to be experimentally informative. In turn An can give some limits to the weak pion coupling 
provided the error can be pushed below 10~ 8 as hoped for the NPDGamma experiment [35(. As for As, being 
sensitive to short range effects, it could carry information on low energy constants of chiral perturbation theory. 
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Appendix A 



The parity conserving and nonconserving effective Hamiltonians (for PNC see [38|]) for the NN and AN vertices 
are given by 



n$% = g P $(f + ^o ,,v qu)T<<.p ll . 



(Al) 



n NNp — V 



,( 2 ) 



hf>T ■ p, + h^pl + ^(3roP° - r • p„) 



7^75-0, 



Hl c Np = i^e-Ths-fT*" ■ {d uPll - dpp,,) + h.c. 



(A2) 
(A3) 
(A4) 



i/PC 
n NNu 



n 



h {1) 

PNC _ '}]<■ 



rCijOlpfiTTjll;, 



(A5) 
(A6) 

(A7) 



Win 



(A8) 



where fp/m p = g*(l + Xp)/2M and is the Rarita-Schwinger vector-spinor field. 



Appendix B 

In this appendix, we give explicitly the reduced matrix elements of the Hamiltonian eq. ([T]) expressed in terms 
of Wigner coefficients. The factor N = */ in the reduced matrix elements (IB1[) - (|B3|) simplifies away in the 
asymmetry observables. The time-reversed adjoint electromagnetic amplitudes differ from the originals only by the 
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phase factor of — (— 1) J , which cancels out in the absolute squares and interference. The reduced NN matrix elements 
for the electric El transition are (with shorthand notations J = 2 J + 1, S = 2S + 1, and L = 2L + 1) 



£:'/(k)=iN(-lf\/3JL'L d 



s' 



V J S' 
1 L d 1 



V 1 L d 




drrU J JZ'(k,r)V K4 (r)6 S 's d . 



(Bl) 



Similarly the reduced NN matrix elements for the magnetic Ml transition are given by 



(-1) 



L' 



3J 



{-lf^L d (L d + l)L d 

M± [(-1)^±(-1) S 



2M 



L' J S' 
1 L d 1 



8s'S d 



S' J V 
1 S d 1 



1 c" 1 

2 D 2 

S d 2 1 



drU J J+\k,r)V Kd (r)6 L , Ld . 



(B2) 



with /i + = /i s and /i_ = /it,. The reduced matrix elements for the magnetic Ml N — > A transitions of particle 1 (for 
isospin transitions 4-> 1 and 1 «-» 2) 



A&(*0= ATV^(-l) 



drU J JZ>{k,r)V Kd (r)6. 



L'+J 



8JS'S d 



I" 


J V 1 


( f 


5" 






s d 1 J 


\ 5 d 







L'L d 



(B3) 



where £ is the degree of freedom for spin: 1/2 for a nucleon and 3/2 for a Delta. In the transitions of particle 2 
cq. (|B3I) has an additional phase of — (— ) Sd+s . Totally the Ml reduced matrix elements are given by -M*f x (k) = 

N£ d J x (k) + A" * x(fc) - For convenience we also define M-% d J \{k) = \M* d J (k). Despite the multiplication factor, the 
general results (|B1I) and (IB2|) reduce to the expressions of ref. [10| for the five lowest amplitudes. 
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